Bistability phenomena in ID polariton wires 
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We investigate the phenomena of the bistability and domain wall propagation in polaritonic 
systems with dissipation provided by the interaction with incoherent phonon bath. The results on 
the temperature dependence of the polariton bistability behavior and polariton neuron switching 
are presented. 



I. INTRODUCTION 

The study of light-matter interactions is interesting 
and important both from the point of view of funda- 
mental fundamental physics and device applications. In 
this context, the structure which attracted a particular 
attention is a semiconductor microcavity (MC) with a 
quantum well (QW) embedded into it. Tuning the en- 
ergy of the excitonic transition in resonance with the en- 
ergy of the photonic cavity mode one can reach a regime 
of strong coupling accompanied by formation of hybrid 
quasiparticles called exciton-polaritons. These half-light, 
half-matter particles exhibit a number of extraordinary 
properties. Due to their extremely small effective mass 
and bosonic nature, polar itons provide an opportunity 
to study various quantum collective phenomena, ranging 
from polariton BECA and Josephson effect to polariton- 
mediated superconductivity^. Strong polariton-polariton 
interactions make it possible to observe several remark- 
able non-linear effects such as superfluidity^, bi- and 
multi-stability, and soliton-like propagation^. 

Although free polaritons are 2D particles, it was re- 
cently noted 8 that it can be interesting to consider ID po- 
lariton wires (see setup in Fig. [1]) provided by the lateral 
confinement of the polaritons in one of the directions^. 
Such wires have the potential to become basic building 
blocks in future spinoptronic devices, including polariton 
Berry phase interferometer— and polariton Datta and 
Das spin transistor 11 . In the nonlinear regime, the bista- 
bility effects in ID polariton channels allow for the real- 
ization of logical circuits based on polariton neuron a 12 i 13 , 
in which a local switching between states propagates 
throughout the wire. 

The above mentioned phenomena require a proper the- 
oretical basis for description of all relevant processes 
which is a non-trivial task. A successful theoretical con- 
sideration should include the effects of interaction with 
a reservoir of acoustic phonons which leads to polari- 
ton thermalization, and the effects of polariton-polariton 
scattering leading to the blueshifts and nonlinear it ies. 
Additionally, it should be taken into account that po- 
laritons have finite lifetime, and the correct description 
of their dynamics should necessarily take into account 
the effects of pump and decay. 

Currently, two main ways have been pursued to de- 
scribe the dynamics of interacting polaritons. Firstly, 
with the assumption of full coherence, mean-field ap- 
proximation gives the Gross-Pitaevskii equation (GPE) 
commonly used for the description of spatially inhomo- 
geneous polariton condensates^. However, while GPE 
includes polariton-polariton scattering it does not de- 



scribe interaction with a phonon reservoir. Assuming 
oppositely that the polaritons are completely incoherent, 
the dynamics in reciprocal space can be described by the 
semiclassical Boltzmann equations^—. Unfortunately, 
this technique fails to describe the real space dynamics 
of inhomogeneous systems. 

Recently, we proposed a novel formalism based on full 
density matrix approach to describe the dynamics of 
the interacting polariton system with dissipation in real 
space and timoiS. The method was applied for the study 
of the propagation of ID polariton droplets. However, 
we neglected the pumping terms and thus consideration 
of the cw regime most relevant from the experimental 
point of view was not done. This paper is devoted to 
the bridging of this evident gap. Here we investigate 
the phenomena of the bistability and domain wall prop- 
agation in polaritonic systems with dissipation provided 
by the interaction with an incoherent phonon bath and 
present the results on the temperature dependence of the 
hysteresis in the polariton system and polariton neuron 
switching. 



II. FORMALISM 

The general formalism for the time evolution of a spa- 
tially inhomogeneous bosonic system without account of 
the coherent pumping terms was developed in our pre- 
vious paper—. We refer to this paper for details and 
derivations and give here only a brief overview, adding 
the detailed consideration of the pumping terms only. 

The state of the whole system is described by the 
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FIG. 1: Schematic representation of polariton wavepacket 
propagation along a ID polariton wire. 
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combined polariton and phonon density matrix p = 
Pph ® ppol (factorization corresponds to Born approxima- 
tion). The phonon part of the system is assumed to be 

time-independent and thermalized, p ph = exp {-/3H ph }, 

while we need to determine the time-dependence of the 
single-particle polariton density matrix in real space. 
The convenient way of its representation is the use of 
the polariton field operators ^(r, t), ^(r, t): 

p(r, r',t) = Tr {ft(r, t)${r', t)p} = (ft (r, t)${r' , t)) (1) 

where the trace is performed by all the degrees of freedom 
of the system. A Fourier transform can be performed to 
work in reciprocal space, which makes the calculations 
easier. 

/ o(k,k / ,t) = (27r) d /^ d J j^-* r,) p(T,T\t)dTdT' = (2) 

where d is the dimensionality of the system (d = 2 for 
non-confined polaritons, d = 1 for the polariton channel), 
L is its linear size, a^, are creation and annihilation 
operators of the polaritons with momentum k. If the time 
dependence of the density matrix in reciprocal space is 
determined, an inverse Fourier transform allows to obtain 
the dynamics in real space straightforwardly. 

The Hamiltonian of the system can be represented as 
a sum of the terms corresponding to various physically 
relevant processes in the system: 



H = Hq + 
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H cp + H ( 
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where 



H = E ^ a l a ^ 



(3) 



(4) 



The effect of the coherent and incoherent parts should 
be treated in different ways. As for the coherent processes 
in the system, their effect can be accounted for using the 
Liouville-von Neumann equation 



i (dtp) 



(co) 



Hco'i P 



(10) 



On the contrary, the incoherent part of the evolution 
is described by the Lindblad equation^ which reads 



(d t pf eco) = - f dt' [H deco (t); [H deco (t'); p(t)}} = 

J — CO 

(11) 

= Sae [2 (H+pH- + H~p H + ) - 

- (H + H~ + H~H+) p-p (H + H~ + H~H + )] 

where the coefficient Sae denotes energy conservation 
and has dimensionality of inverse energy and is in the 
calculation taken to be equal to the broadening of the po- 
lariton state^i. The terms H + and H~ correspond to the 
processes when the thermal reservoir particle (phonon or 
other reservoir boson, see below) is created or destroyed. 

The effects of the terms corresponding to the polariton- 
phonon and polariton-polariton interactions was in de- 
tail considered in Ref. 19. The corresponding equations 
for the elements of the density matrix and their deriva- 
tion can be found there, and we do not reproduce them 
in current manuscript. For polariton- polariton scat- 
tering, these equations reproduce an analog of Gross- 
Pitaevskii equation written for the density matrix, and 
for polariton- phonon interaction they are generalizations 
of semiclassical Boltzmann equations. 

Here we only consider the terms provided by coherent 
and incoherent pumps, which should be added to the 
dynamic equations of Ref. [l9 to make our consideration 
complete. 



corresponds to free polariton propagation, 
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to polariton-polariton scattering, 



(5) 



H P h = £(q)4 +q a k (& q + &- q )> (6) 

k,q 

H cp to coherent laser pumping, Hdcp to incoherent pump- 
ing and finite polariton lifetime (see expressions for these 
two terms below). 

In the above formulae defines the dispersion of the 
free polaritons and the quantities U and D correspond 
to the polariton- polariton and polariton- phonon scat- 
tering. 

The Hamiltonian ([3]) can be separated into the sum of 
the coherent part and the part introducing decoherence, 



H = Hn 



decoi 



io — Hq - 
H-deco 



■ Hpol 

= H p h 



dcp 



(7) 
(8) 
(9) 



A. Coherent pumping 

We start from the case of the coherent pump. Its phys- 
ical meaning is coupling to an electric field with a well 
defined phase, provided e.g. by an external laser beam. 
Mathematically, the corresponding Hamiltonian can be 
introduced as 



(12) 



The coefficients pk are Fourier transforms of the pumping 
amplitude in a real space by p(x, t), which can be cast as 



p(x,t) = P(x)e^ x e 



(13) 



where P(x) is the pumping spot profile in real space, k p is 
an in-plane pumping vector resulting from the inclination 
of the laser beam as respect to the vertical and uo p is the 
pumping frequency of the single- mode laser. 

Let us now check the effect of H cp on the evolution 
of the polariton density matrix. Insertion of the Hamil- 
tonian (fT2|) into Liouville- von Neumann equation yields 
the following result: 
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ft(4flk) = -^lmK(flk)} (14) 

^(o^ Ok 2 > = ^(4K> -Pk 2 (a kl )*) 

One sees, that the equations for matrix elements con- 
tain a new quantity- an average value of the annihilation 
operator of the polariton field, which is nothing but the 
order parameter (also called a macroscopic wavefunction) 
of the system. One thus needs to obtain the dexpression 
for this quantity to close the system of the equations. 
Straightforward derivation gives: 



<9t(a k ) = --Pk - ^^k(«k) (15) 
~ \ U Yl p ( k<2 ' k 2 - P) ( a k+ P ) 

k2,P 

+ ( E W(q)(p(k + q,k + q)-<) 

+ E W ^)(-P( k + q, k + q) - < - 1) ) (a k > 

where U is a matrix element of polariton- polairton inter- 
action and W(q) are transition rates of phonon- assisted 
processes (See^ for the definition of the parameters). 



B. Incoherent pumping and lifetime 

By incoherent pump we mean the exchange of parti- 
cles between the polariton system and some incoherent 
bosonic reservoir, whose nature depends on the pumping 
scheme. Usually, this will be an ensemble of incoher- 
ent excitons created either by an electrical pump or by 
incoherent optical excitation or a reservoir of the exter- 
nal photonic modes providing the leakage of the photons 
from the cavity. The corresponding Hamiltonian written 
in Dirac representation reads 

H dcp = J2 K ^ kV (^-^W^&k/ + H.c. = (16) 

k,k' 

_ rx+ _i_ zj 

icp * icp 

where 6k is a secondary quantization operator corre- 
sponding to the bosonic reservoir in question, if (k, k') 
are constants characterizing the coupling between the po- 
lariton system and the reservoir The introduction of this 
Hamiltonian into Lindblad equation leads to the stan- 
dard terms, whose derivation can be found elsewher e 20 ! 23 : 



dt(<a k /) = ik<W - 7^ (7k + 7k')( a ]k a k') (17) 

where the terms 2k and 7k denote the intensity of the 
incoherent pump and broadening of the polaritonic lev- 
els connected with lifetimes of the polariton states, 7k = 



Hr^ . They are usually taken as phenomenological pa- 
rameters, but can be connected with the quantities en- 
tering in Hamiltonian ([T7|) : 

/k = \ E l^(k,k')| 2 ^(k) _ E'(k'))n w , (18) 

k' 

7k = ^|^(k,k')| 2 Wk)- J E'(k')) (19) 

k' 

where are the occupancies of the bosonic reservoir. 

In our further consideration we will consider only the 
case of the coherent pump, thus putting all 1^ = and 
retaining in the resulting equations only the terms corre- 
sponding to the lifetime, which corresponds to the case 
of the empty bosonic reservoir , = for all the k. 



III. RESULTS AND DISCUSSION 

To get a full system of equations for the dynamics of 
the polariton system with pump and decay, one should 
combine the equations derived in Ref. Qj)] with expres- 
sions ([H |h([TI^ . (fTT]) . The resulting formalism is suitable 
both for describing 2D polaritons and polaritons confined 
in ID channels. The consideration of the former case, 
however, needs powerful computing facilities and in the 
present paper we focus only on consideration of the latter 
one. 

We considered a 2 jim wide polariton channel in a mi- 
crocavity with an active region based on InAlGaAs alloys 
with Rabi splitting 15 meV. The matrix elements of the 
polariton- polariton and polariton- phonon interactions 
were estimated using the standard formula o 24 i 25 . 

The first phenomenon we modeled was the effect of 
the thermalization in the polariton system provided by 
polariton- phonon interaction on its bistable behavior. 
It is well known that if the polariton ensemble is fully 
coherent and its dynamics is described by the Gross- 
Pitaevskii equation containing coherent pumping and 
lifetime terms, and if the energy of the pumping laser lies 
slightly above the bottom of the lower polariton branch, 
the dependence of the concentration of the polaritons on 
pump intensity is described by an S-shaped curve, char- 
acteristic for systems revealing the effects of bistability 
and hysteresis 2 ^. Such a behavior is due to the polariton- 
polariton interactions which introduce nonlinearity into 
the system. On the other hand, in the approach based on 
the semiclassical Boltzmann equations corresponding to 
the limit of strong decoherence the bistability is absent 
and the dependence of the occupancy of the ground state 
on the pump intensity is described by a single defined 
threshold function 1 -. One can expect that a transition 
between these two regimes should occur if one raises the 
temperature in the system which leads to the intensifi- 
cation of the polariton- phonon interactions and deco- 
herence in the system. This was indeed observed in our 
calculations. 

The computational results are presented in Fig. [2j The 
system is pumped with a spatially homogenious laser 
beam oriented perpendicular to the QW (i.e. at k = 
in /c-space) at a slightly higher frequency than the k = 
polariton frequency (0.24 meV detuning). At low pump 
intensity, the pump frequency is not in resonance with 
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FIG. 2: The dependence of the polariton concentration on 
spatially homogeneous cw pumping intensity for various tem- 
peratures. At 20K, the hysteresis curve is quite wide (black). 
At 30K, the hysteresis curve is much narrower, while at 50K 
it has disappeared completely. 



the condensate. Consequently, the condensate occupa- 
tion remains fairly low. As the pumping is increased, the 
polariton energy is blueshifted into resonance with the 
pump and there is a sudden jump in the occupation of 
the k = polariton state at some characteristic pump in- 
tensity Iq. If one then decreases the pump, the polariton 
occupancy jumps down at different value I\ < Iq which 
corresponds to the hysteresis behavior. However, the in- 
crease of the temperature leads to the intensification of 
the phonon scattering, and bistable behavior becomes 
less and less pronounced: the hysteresis area narrows and 
is quenched completely above some critical temperature 
T c & 4QK. This correspond to the transition between the 
Gross-Pitaevskii and Boltzmann regimes in the polariton 
system. A similar phenomenon was earlier predicted for 
microcavity- based terahertz emitting device^. 

The effect of bistability can form a basis for cre- 
ation of a variety of devices based on nonlinear polari- 
ton transport^. Among them are polariton neuronsi^- 
the building blocks of polariton based optical integrated 
circuits^ utilizing the phenomenon of domain wall prop- 
agation in bistable systems. The underlying idea is fol- 
lowing. Imagine that the polariton system is driven by a 
spatially homogeneous background cw pump with inten- 
sity corresponding to the bistable regime and the steady 
state occupancy of the system corresponds to the lower 
branch of the S- shaped curve. Then a short localized 
pulse is applied in the middle of the wire. Its intensity 
should be enough to send the condensate locally to the 
upper branch of the S-shaped curve. Due to diffusion, 
the polariton wavepacket spreads to the neighboring re- 
gions and switches them to the upper branch. This way, 
the area of high occupancy steadily expands. This phe- 
nomenon is analogous to the propagation of the domain 
wall in ferromagnetic materials. It should be noted that 
although the polar itons have finite lifetime, this does not 
limit the length of signal propagation in a polariton neu- 
ron, and the signal keeps propagating as long as the back- 
ground cw pumping persists. 

As bistability in the polaritonic system strongly de- 
pends on temperature as was discussed above, the same 
should be true for the domain wall propagation in the 
polariton neurons. This can indeed be seen from Fig. [3] 
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FIG. 3: Neuron behavior and domain wall propagation for 
(top to bottom) T=15K, 20K and 25K. The plots show the 
polariton concentration in real space (horizontal axis) and 
time (vertical axis). The system is pumped non-resonantly 
by spatially homogenious cw laser, at 70-80ps a switching 
pulse arrives. At 15 K one sees the propagation of the dis- 
tinctive domain wall which becomes more smeared and fully 
disappears at 25 K due to that there is no bistability at this 
pump intensity and temperature. 



representing the results of our calculations. At 15K, the 
switching is clear and one can easily detect the propa- 
gation of a sharp well defined domain wall. At 20K and 
using the same pump intensity as before, the switching 
is less pronounced, partly because the higher intensity 
is lowered and partly because the spreading of the high 
intensity is slower, which leads to the washing out of the 
domain wall. At 25K there is no bistability for this pump 
intensity anymore, and the system goes back to its low 
population state after the application of the pulse. 

Finally, the effect of phonon scattering on pure de- 
phasing in the system was considered. We pumped the 
system with a coherent pulse having Gaussian profile in 
the real space. For various temperatures of the system 
state we investigated in the steady state the spatial pro- 
files of both total polariton density and its coherent part 
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FIG. 4: Coherent fraction at (top to bottom) T=10K, 30K 
and 60K. Pumping at k — 0, no detuning. The solid lines show 
the total population, the dashed lines the coherent population 
and the dotted ones the incoherent population. 



determined as \i(j(x)\ 2 , where 

/+oo 
(a k )e ikx dk 
-oo 

The results are shown in Fig. |4j One can see that at low 
temperatures the coherent fraction is quite large at the 
center of the coherent pumping spot, but it dramatically 
decays outside it. However, as temperature is increased, 
the density profile gets more spread out over the wire, 
forming an almost constant background density made up 
of the decoherent population. 



IV. CONCLUSION 

In conclusion, we have considered the effects of coher- 
ent pumping and finite lifetime in a polaritonic system 
accounting for all physically relevant processes in the sys- 
tem. We applied our theory for the consideration of non- 
linear polar iton propagation in a ID polar iton wire. We 
have shown that the increase of temperature dramati- 
cally affects such processes as bistability switching and 
domain wall propagation in polar iton neurons. 
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